Parton rescattering effect on particle production in ultra-relativistic p+p collisions 
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The parton rescattering effect on the charged particle production in ultra-relativistic p+p colli- 
sions is studied by the parton and hadron cascade model, PACIAE, based on PYTHIA. We have 
calculated charged particle pseudorapidity density (dN c h/dri) at mid-rapidity and the pseudora- 
pidity distribution in inelastic and non-single diffractive p+p collisions at ^=200, 900, 5500, and 
14000 GeV with the PYTHIA and PACIAE models. The calculated results of ^=900 GeV are well 
compared with the ALICE data. The calculated dN c h/drj as a function of center-of-mass energy 
well meets with the experimental data as well. Comparing the PYTHIA results (without parton 
rescattering) with the PACIAE results (with parton rescattering), it turned out that the parton 
rescattering effect plays an important role and this effect increases with increasing center-of-mass 
energy. 

PACS numbers: 25.75.Dw, 24.10.Lx 



I. INTRODUCTION 



The jets, as remnants of hard-scattered quarks and 
gluons, have been investigated to explore the proper- 
ties of partons Studying the hot QCD matter via 
jet quenching (partons lose energy on their way through 
the QCD medium) has now become an important role 
of jets and has been realized at RHIC 0, Q| . Experi- 
mentally extracting the partonic observables from final 
state hadrons with reconstruction method is now highly 
interesting [4|, |5j. 

Another way to study the properties of partons is 
introduced recently. The pr distribution of effective 
u (d) and s quarks are extracted from the ratio of 
E(p T /3)/<j}(p T /2) and fi(pr/3)/0(pr/2) in the Au+Au 
collisions at v /j/v r /v=200 GeV , respectively [1,0] ■ In this 
extraction it is assumed that the hadron's pt is com- 
posed of its effective constituent quark's pr/n (n is the 
number of effective constituent quarks in the hadron) 
and the different quark and anti-quark have the same 
Pt distribution. 

The measurement of elliptic flow parameter V2 relies 
on the analysis techniques which require high event mul- 
tiplicity. Therefore, no report about the V2 measurement 
in p+p collisions is published by now because of low mul- 
tiplicity. However, in recent papers 0,0, it is argued 
that the elliptic flow parameter in high multiplicity p+p 
collisions at the LHC energy (y / 5=14000 GeV) may be 
measurable. Similarly, one could expect the formation 
probability of quark-gluon mater (QGM) in the p+p col- 
lisions at RHIC energy (^5=200 GeV) is very low [lj. 
The parton rescattering effect on the final state hadrons 
in p+p collisions at RHIC energy is weak as well [121 ]. 
However, those probability and effect may be visible in 
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the p+p collisions at LHC or higher energy [T3 |. 

The study of ultra-relativistic p+p collisions is a 
door sill toward the nucleus-nucleus collisions at the 
same energy. Because the nucleus-nucleus collision can 
be decomposed into nucleon-nucleon collisions and the 
nucleon-nucleon collision can be well described by the 
perturbative QCD (pQCD) due to its relatively clear 
and simple physics. The results of p+p collisions are 
a worthy reference for nucleus- nucleus collisions. Thus 
the ALICE first successful measurements for p+p col- 
lisions at -^=900 GeV [l3[ not only demonstrate "the 
LHC and its experiments have finally entered the phase 
of physics exploitation" but also light the way studying 
Pb+Pb collisions at LHC energies. It also means the 
copious measurements in p+p collisions at LHC energies 
are forthcoming and will supply opportunities for the 
investigations of the issues mentioned above. 

As a dedication to the ALICE first measurements, in 
this paper the parton and hadron cascade model PA- 
CIAE, based on PYTHIA, is used to analyze the AL- 
ICE first measurements, especially to explore the parton 
rescattering effect on the final state hadron in the ultra- 
relativistic p+p collisions. This investigation is not only 
a quick response to the ALICE first measurements but 
also a well preparation for the further studies in p+p 
and Pb+Pb collisions at LHC energies. 



II. MODELS 

The parton and hadron cascade model, PACIAE 
is based on PYTHIA [3. PYTHIA is a model for high 
energy hadron-hadron (hh) collisions. In the PYTHIA 
model a hh collision is decomposed into the parton- 
parton collisions. A hard parton-parton collision is de- 
scribed by the lowest leading order perturbative QCD 
(LO-pQCD) parton-parton interactions modified by the 
parton distribution function in a hadron. The soft 
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parton-parton collision is considered empirically. Be- 
cause the initial- and final-state QCD radiations are con- 
sidered in parton-parton interactions, the consequence of 
a hh collision is a partonic multijet configuration com- 
posed of di-quarks (anti-diquarks), quarks (anti-quarks), 
and gluons. This is followed by the string construction 
and the string fragmentation. So one obtains a final 
hadronic state for a hh (p+p) collision. 

The differences between the PACIAE model and 
PYTHIA model, for p+p collision, are as follows: 

• The string fragmentation mentioned above is 
switch-off temporarily in the PACIAE model. 
Therefore, if the di-quarks (anti-diquarks) are bro- 
ken randomly into quarks (anti-quarks) the conse- 
quence of a hh (p+p) collision is a configuration 
of quarks, anti-quarks, and gluons. As such the 
partonic initialization stage for a p+p collision is 
realized. 

• Then the rescattering among partons is introduced 
and performed by the Monte Carlo method until 
partonic freeze-out (no more parton-parton inter- 
action at all). In this parton rescattering stage the 
LO-pQCD parton-parton interaction cross sections 
(2 — > 2) are employed [l6l |. 

• In the next hadronization stage, the partonic mat- 



ter formed after parton rescattering is hadronized 
by the Lund string fragmentation regime [TEj 
and/or coalescence model [3~3 | . 

• The consequent hadronic matter suffers hadronic 
rescattering. We deal with hadronic rescattering as 
usual two-body collision [l7| , until hadronic freeze- 
out (hh collision pair is exhausted). This is the 
hadronic rescattering stage. 



III. PARTON RESCATTERING EFFECT ON 
PARTICLE PRODUCTION 

As we aim at the physics behind the experimental data 
rather than reproducing the data, in all calculations the 
model parameters are fixed at their default values, ex- 
cept the k factor, considering the higher order and un- 
perturbative QCD corrections, is assumed to be equal to 
3. Both inelastic (INEL) p+p and non-single diffractive 
(NSD) p+p collisions are calculated with and without 
parton rescattering at ^5=200, 900, 5500, and 14000 
GeV. The calculated results without parton rescattering 
will be indicated as PYTHIA and the results with parton 
rescattering as PACIAE. In addition, a similar calcula- 
tions are also performed for the p+p collisions at s/s= 
900 GeV. 



TABLE I: Charged particle pseudorapidity density at mid-rapidity (|?7| < 0.5) calculated by the PYTHIA and PACIAE models 
for p+p collisions at various center-of-mass energies. The ALICE p+p and UA5 p+p data [l|, Ql| at y^—QOO GeV are given 
as well. In the ALICE data the first error is statistical error and the second is systematic one. 

P+P P+P 



ALICE 320 tune 1 PYTHIA PACIAE UA5 PYTHIA PACIAE 

INEL 3.10 ± 0.13 ± 0.22 2.46 2.59 2.80 3.09 ± 0.05 2.70 2.94 

NSD 3.51 ± 0.15 ± 0.25 3.02 3.15 3.33 3.43 ± 0.05 3.30 3.50 
1 The PYTHIA results of the Perugia (320) tune taken from [21]. 



Table I gives the calculated charged particle pseudo- 
rapidity density at mid-rapidity (h) | <0.5), dN c h/dr]. to- 
gether with the ALICE data [lj| and UA5 data [II]. 
One sees in this table that: 

• The PYTHIA results of dN ch /dr) in this paper are 
within the range of three PYTHIA tunes [13| : from 
2.33 01 to 2.99 [H for INEL and from 2.83 [H| to 
3.68 H for NSD. Our PYTHIA results are close 
to the Perugia (320) tune [Hj]. 

• Both the ALICE and UA5 data are better repro- 
duced by PACIAE than PYTHIA. This means the 
parton rescattering effect is important in reproduc- 



ing the experimental data. 

• The difference in charged particle pseudorapidity 
density between NSD p+p and p+p collisions at 
Vs=900 GeV is about 4.5% in the PYTHIA calcu- 
lations and 4.9% in PACIAE. They are larger than 
the range of 0.1-0.2% resulted from extrapolating 
the CERN ISR measurement [H to yfi =900 GeV. 

Figure [1] gives ALICE data of charged particle pseudo- 
rapidity distributions in p+p collisions at ^fs =900 GeV 
[H| together with the results from PYTHIA and PA- 
CIAE calculations. We see in this figure that the ALICE 
data are better reproduced by PACIAE than PYTHIA. 
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FIG. 1: (Color online) Charged particle pseudorapidity dis- 
tributions in INEL and NSD p+p collisions at ^=900 GeV. 
The full squares and triangles are ALICE INEL and NSD 
data, respectively, taken from [l3| . The open squares and 
triangles are the PACIAE results. The dashed and dotted 
curves are the PYTHIA results. 



Both PYTHIA and PACIAE results show the similar 
shallow valley at r\ ~0 to the ALICE data, and the bet- 
ter symmetry relative to rj=0. 

The calculated center-of-mass energy dependence of 
the charged particle pseudorapidity density at mid- 
rapidity is compared with the ALICE data in Fig. [21 
One sees in this figure that the power-law fit is bet- 
ter reproduced by PACIAE rather than PYTHIA both 
for INEL and NSD. The PYTHIA results are below the 
corresponding power-law fit lines for both INEL and 
NSD. The PACIAE NSD results deviate from the cor- 
responding power-law fit lines at y/s ^2500 GeV. Then 
this deviation increases with increasing center-of-mass 
energy. The role of parton rescattering also increases 
with increasing center-of-mass energy. For instance, the 
difference between PACIAE and PYTHIA NSD results 
(in percentage) increases from 5.1% at -^=200 GeV to 
15.4% at 5500 GeV and then to 24.2% at 14000 GeV. 

We compare the charged particle pseudorapidity dis- 
tributions calculated by the PYTHIA and PACIAE mod- 
els for NSD p+p collisions at ^=200, 900, 5500, and 
14000 GeV in Fig. [3] In this figure the open and full 
symbols are the PYTHIA and PACIAE results, respec- 
tively. The circles, squares, triangles, and diamonds in- 
dicate the results in ^7=200, 900, 5500, and 14000 GeV 
p+p collisions, respectively. One sees the following fea- 
tures in this figure: 

• The charged particle pseudorapidity density at 
mid-rapidity and the width of center rapidity 
plateau increase with increasing center-of-mass en- 
ergy monotonously. 



FIG. 2: (Color online) Charged particle pseudorapidity den- 
sity at mid-rapidity (\rj\ <0.5) in p+p collisions as a func- 
tion of the center-of-mass energy. The triangles and circles 
are the PYTHIA and PACIAE results, respectively. The full 
and open symbols indicate the NSD and INEL p+p collisions. 
The dashed and solid lines indicate the fits to experimental 
data (y/s < 1900 GeV) by a power-law dependence on the 
energy for INEL and NSD interactions, respectively [l3| ]. 



• In the center pseudorapidity region the PACIAE 
results are larger than PYTHIA. However, one sees 
the opposed situations in the outer pseudorapidity 
region. This trend increases with increasing center- 
of-mass energy. That means the parton rescatter- 
ing effect (the QGM formation probability) in p+p 
collisions increases with increasing center-of-mass 
energy as well. 



IV. CONCLUSIONS 

Using the PYTHIA model (without parton rescat- 
tering) and the parton and hadron cascade model PA- 
CIAE (with parton rescattering) we have investigated 
the charged particle productions in INEL and NSD p+p 
collisions at Vs=200, 900, 5500, and 14000 GeV as well 
as p+p collisions at \/s—900 GeV. 

The calculated charged particle pseudorapidity den- 
sity at mid-rapidity (\rj\ <0.5) and the pseudorapidity 
distributions in p+p and p+p collisions at y / s=900 GeV 
well reproduce the corresponding ALICE and UA5 data 
0i EH- The calculated dN c h/drj at mid-rapidity as a 
function of center-of-mass energy well meets with the 
power-law dependence lines fitted to the experimental 
data in the yfs <1900 GeV region (l3j |. Above this region 
the deviation of PACIAE NSD results from NSD power- 
law dependence line increases with increasing center-of- 
mass energy. The calculated charged particle pseudo- 
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rapidity distributions, similar to the ALICE data, show 
a nearly plateau structure at mid-rapidity with a shal- 
low valley at center pseudorapidity. That plateau width 
increases with increasing center-of-mass energy. 

The parton rescattering effect plays an important 
role in reproducing the experimental data. The par- 
ton rescattering effect becomes larger and larger with 
increasing center-of-mass energy. This also means the 
QGM formation probability in the p+p collisions in- 
creases with increasing center-of-mass energy. Therefore 
the QGM is possible to be formed in the early stage of 
p+p collisions at ALICE or higher energy. 



FIG. 3: (Color online) Charged particle pseudorapidity dis- 
tributions in NSD p+p collisions at y/s—2QQ (circles), 900 
(squares), 5500 (triangles), and 14000 GeV (diamonds) cal- 
culated by PYTHIA (open symbols) and PACIAE (full sym- 
bols) models, respectively. 
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